Study of Multiple Rounds of Error Correction in Solid State NMR QIP 
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Methods to control errors will be essential for quantum information processing. It is widely 
believed that fault-tolerant quantum error correction is the leading contender to achieve this goal. 
Although the theory of fault-tolerant quantum error correction is very well understood, experimental 
implementation has been lagging. We study the feasibility of implementing repeated rounds of 
quantum error correction with refreshed ancillas in solid state nuclear magnetic resonance (NMR). 
In particular we study the procedure proposed for extracting entropy that consists of two stages; 
an error correcting code optimized to function at finite temperature, and an implementation of 
heat-bath algorithmic cooling to refresh the ancilla qubits. Two algorithms are presented which 
implement this method, one for performing tests on 4 qubits, the other for practical implementation 
on 6 qubits. The effects of imperfect implementation are examined in both the error correction and 
refreshing stages. 



I. INTRODUCTION 

It is possible to encode information onto systems which 
are governed by the laws of quantum mechanics. The ex- 
tra computational freedom granted by such an approach 
can enable exponentially faster algorithms for a growing 
set of tasks. Included among these are the factoring of 
large numbers [1 , database search [2], and simulation of 
physical systems governed by quantum mechanics [3J. 

A challenge to the implementation of quantum algo- 
rithms in physical systems, with which this manuscript 
is concerned, is the process of decoherence. Quantum 
states are inherently delicate; processes beyond experi- 
mental control drive quantum systems toward classical 
behaviour. This presents a challenge to the implementa- 
tion of quantum computing, whose efficacy depends on 
the presence of non-classical behaviour. 

There are two principal means of reducing or eliminat- 
ing these errors. The first is to further isolate the quan- 
tum computer from the external environment, reducing 
the unwanted interactions directly. This is the tactic 
employed by topological quantum computing [53], which 
uses topologically-protected states to store information. 
The second method, discussed in this manuscript, is the 
use of error correction protocols (as seen in Figure [lj . 
These are based on encoding a one-qubit state into mul- 
tiple qubits, allowing an error process to act upon the en- 
coded state, and decoding/correcting the original state. 

Some elements of fault-tolerant quantum error correc- 
tion have already been implemented. In many systems, 
phase errors are dominant and thus the 3 qubit phase 
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correcting code ([H], [H], [2D], EU) and a 5 " bit code 
[22] have been implemented to reduce the noise. 

In the case of NMR, it is not possible to make projec- 
tive measurement adding difficulty to the task of refresh- 
ing qubits for more rounds of error correction. It is how- 
ever possible to use algorithmic cooling as suggested by 
([15] , [16]) and first demonstrated in [7]. More recently, 
another building block for fault tolerance was demon- 
strated: magic state distillation (sec 24J for theoretical 
background, |25| for details of implementation). 

In order to practically implement quantum computing, 
it is necessary to have the ability to perform multiple iter- 
ations of error correction, so that states can be preserved 
against arbitrary noise over arbitrary time scales. This 
arises from the presence of second-order error terms. 

A common assumption underlying the design of error 
correcting codes is that ancilla qubits (qubits employed 
in encoding in addition to the qubit on which the message 
is stored) are highly polarized, typically initialized in the 
state |0). These ancilla qubits, in practice, are in short 
supply and need to be 'refreshed' to high polarization in 
order to run multiple iterations of error correction. 

In section II of this manuscript, the proposed experi- 
mental setup is discussed, in order to provide context to 
the algorithm detailed in the following sections. In sec- 
tion III of this manuscript, a modified version of the tra- 
ditional three-bit error correcting code is presented, with 
the experimental implementation of multiple rounds of 
error correction as the eventual goal. In Section IV, a 
known mechanism for transferring entropy is analysed, 
and the effects of imperfect control in experimental error 
correction are discussed. 
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FIG. 1. Sketch of a generic error correction procedure, decom- 
posed into three phases. First, the one-qubit basis states (typ- 
ically labelled |0) and |1)) are encoded (E) into two distinct 
n-qubit states. Then the encoded state is transmitted (T) 
through a noisy channel, or stored in an unreliable medium. 
In the third phase, the state is decoded (D), providing an 
output state with improved fidelity to the original state. The 
ancilla qubits, incorporated in the encoding phase, are then 
disposed of. 



II. EXPERIMENTAL SETUP 

In order to provide context for the code presented, we 
first detail the experimental situation for which this code 
is designed. This does not preclude the code presented 
in later sections from being used in other quantum com- 
puting architectures, provided that similar resources are 
present. 

Single-crystal solid-state NMR employs a large ensem- 
ble of identical molecules with fixed orientation on a crys- 
tal lattice. Each nucleus in a molecule precesses about 
an axis defined by the large external magnetic field at its 
Larmor frequency w, which is the product of a 'chemi- 
cal shift' term (1 + S) which arises from local electrons 
screening the external magnetic field, the gyromagnetic 
ratio of the nucleus in question 7, and the magnetic field 
Bq . The state of the spins is measured by analysing the 
response of the system to an rf pulse. Two-qubit opera- 
tions are obtained through dipole-dipole coupling. These, 
together, allow for universal control, and the large num- 
ber of identical quantum information processors produces 
a detectable signal. The advantages of this implementa- 
tion of quantum computing are described further in [¥], 

The physical system employed in [5], [5] and [7] mer- 
its discussion, since it is similar in many respects to 
the system envisioned for the experimental implementa- 
tion of the procedures discussed below. The architecture 
used is single-crystal solid-state nuclear magnetic reso- 
nance in which malonic acid (the structure of which is 
shown in Figure [2| with three 13 C is present in the crys- 
tal sample at ~ 3% abundance. The remainder of the 
molecules have natural abundance 13 C. These unlabcllcd 
molecules prvent unwanted interactions between labelled 
molecules by ensuring that the probability of 2 triply- 
labelled molecules neighbouring one another in the crys- 
tal lattice remains negligibly small. The algorithms pre- 
sented in the sections below will require 4 and 6 qubits, 
and a high-polarization bath. Larger molecules will be 
required in order to perform these experiments. 

In NMR quantum computing, the dominant source of 
error is 'dephasing'. Any process that involves uncon- 
trolled rotations about the quantization axis (the z-axis) 
results in a dephasing error. In NMR, this uncontrolled 




FIG. 2. Malonic acid. This molecule encodes 3 qubits in 
NMR experiments, and is the molecule used in [7]. A malonyl 
radical can be formed by removing one of the protons; the 
remaining electron can be coupled to the quantum system of 
interest. 



rotation is either the result of a difference in the mag- 
netic field to which each individual spin is subjected, or 
the result of parasitic coupling between qubits and non- 
qubit particles present in the crystal. We represent this 
as a quantum channel Ad op hasing involving the random 
application of the Pauli Z operator with probability p. 
This channel can be expressed in the operator-sum rep- 
resentation: 



A 



dephasing 



1 



-pi, Vp^} 



A( P ) = ^A jP A] 



= (! -p)p + pZpZ 



(1) 



(2) 



Here, I is the identity operator. We can correct such 
errors with a variation on the traditional 3-bit error cor- 
recting code ([8], [9]), that translates dephasing errors 
into 'bit-flip' errors, shown in Figure [3] This error cor- 
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FIG. 3. The 3-bit error-correcting code, modified to pro- 
tect against dephasing errors. The first qubit is prepared 
in an arbitrary state \ip) = a |0) +/3|1). The other two 
qubits are prepared in the state |0). The state on the first 
qubit is 'encoded' to a 000} + /3|111). Then, each qubit 
is subjected to the Kraus map defined in Equation [I] The 
most likely outcomes (7" 183 has no effect, or flips a single 
qubit) map the 3-qubit codeword into two orthogonal sub- 
spaces; 000) -> {j000}, |001), jOlO), 1 100}} and |111} -> 
{jlll} , 1 110} , |101> , |011}}. These errors can then be dis- 
tinguished and corrected, to yield an error probability that 
scales as p 2 instead of p, since two independent errors must 
occur in order to induce an error in the decoded state. 

recting code employs Hadamard gates to transform the 
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dephasing channel into the more familiar bit-flip channel. 
The Hadamard gate transforms the states {|0) , |1)} to 
equal superpositions of |0) and |1), its matrix form is 
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In order to determine the effect of the process 
H Adephasing-ff , we note two properties of the Hadamard 
operator: 



H 2 = 1: HZH = X 



(4) 



HA. 



,1. phasing-^ = y/pXj 

= Abit-flip- (5) 



Another important feature of NMR quantum comput- 
ing is the ability to form radicals (molecules with one pro- 
ton missing) with electron spins that can be controlled 
and coupled to the nuclear system through ENDOR [ID] . 
If electron-nucleus interactions can be controlled, one or 
more electrons in each molecule can be used as a heat 
bath. Even though all components of a given molecule 
are at the same temperature, the polarization in the spin 
degree of freedom of particles inside a molecule may dif- 
fer, since the polarization depends also on the magnetic 
field, and most importantly on the gyromagnetic ratio 7 
of the particle in question: 



tanh 



(6) 



Since the gyromagnetic ratio of the electron is ~ 10 3 
times greater than the gyromagnetic ratio of 13 C, the 
electron is more polarized than any nucleus at thermal 
equilibrium. A transfer of polarization from electrons 
to nuclei is, effectively, an exchange of entropy or 'heat' 
between the two particles. This capacity for entropy ex- 
change will be used to polarize 13 C beyond the bath po- 
larization. At temperatures on the order of 1 K, the 
electron polarization is of order 1, allowing for unit po- 
larization on the nucleus as well. 



III. OPTIMAL ERROR CORRECTION 

The existing three-qubit code has been developed with 
the underlying assumption of polarized ancillary qubits 
which can be quickly exchanged. This code is presented 
in Figure [3] For comparison, the error channels whose 
effects must be corrected are presented in Figure [2] Note 
that, in addition to the channel which is corrected by the 
traditional code, there is a second channel which maps 
the ideal (polarized) ancilla qubits to thermal states. In 
this section of the manuscript, we review the methods 
used to optimize error- correcting codes for experimental 
implementation. This section concludes with the pre- 
sentation of an improved error correcting code analysed 



throughout the remainder of the manuscript and limits 
on the degree of mixture that is tolerable for such a code 
to be useful. 



A. Background 

In order to optimize error-correcting codes, we require 
an objective function to quantify their performance; a 
measure of the similarity between the input and output 
states. The 'channel fidelity' Fc, a special case of Schu- 
macher's entanglement fidelity [TT] suits this purpose: 

F C (A) = (0| (A® 1) [|fi) |fi> = ~£ |tr(A 3 )| 2 ; (7) 

3 

where |fi) = — = ^ \jj) 

Here, A is a quantum channel representing all of the 
error-correction operations (encoding, error, decoding), 
with Kraus operators A,-. |f2) is an entangled state with 
2n qubits, where the error correction is to take an n-qubit 
state as the input. By extending the quantum system, 
the channel fidelity provides a concise means of averaging 
fidelity over all n-qubit states. 

With the metric in Equation [7] we can now define 
a criterion for useful error correction. We require that 
F C {T> o T o £ o M o _4)(where V oT ofoMo^isa 
composite process consisting of appending ancilla qubits 
in the state |00) (00|, A[p] = (1® |00) (00|)p(l® jOO) (00|), 
mixing the states on the ancilla (A4 , see Equation [TT] ) , 
encoding the resulting 3-qubit state (£), transmission (/) 
and decoding (T>)) is greater than Fc{T). Fc{T) is triv- 
ial to evaluate: 
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The channel fidelity also provides a means of optimizing 
the encoding and decoding operations shown in Figure 
[TJ based on 'semi-definite programming' (SDP), a well- 
studied class of optimization problems |12j . 

The channel fidelity of the composite channel (referred 
to as r, see Equation [7]) can be written as a vector inner 
product by stacking the columns of the matrices T> ® 
1[|0) and (To£oM oA)„ ® where A* is 

the Schrodinger dual of a channel (A* [p] = J^j A] pAj , 

cf. Equation B. The CPTP constraint Efc ^I^fe = 1 
can be expressed in SDP. If we do not apply additional 
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constraints, the optimal encoder found by this procedure 
incorporates the following CPTP channel: 
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This channel, when implemented before the encoding 
unitary, returns the ancilla to the state 1 00) with 100% 
fidelity. Physically, this is equivalent to allowing the an- 
cilla to come to equilibrium with a heat bath at zero 
temperature |13) . and it cannot be implemented using a 
unitary operator. In order to find an optimal unitary, we 
further restrict the SDP to unital channels (those that 
obey the constraint J^k^^k = !)• Upon finding the 
channel fidelity of the optimal channel subject to these 
constraints, a unitary operator which replicates this fi- 
delity is found by inspection. According to the scheme 
in [T3], SDP can be iterated to find a code which has 
a globally optimized encoder, and a globally optimized 
decoder, and is, on the whole, at least locally optimal. 



B. Model Channel 

In contrast with the ideal situation typically studied, 
each of the ancilla qubits in an NMR experiment is pre- 
pared in a mixed state. This is equivalent to the action of 
a second error channel which occurs before the encoding 
operation. A mixing channel which works on all inputs 
is identical to the depolarization channel [13] ; a simpler 
channel can be used when the input state is |0): 



M = yi-|i; M 1 = ^|a ; M = {M 0) M 1 }. 
Treating the presence of entropy in the ancilla qubits 
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FIG. 4. A generic 3-qubit error correction procedure, with 
the encoding and decoding procedures replaced by unknown 
processes. The presence of initially mixed ancilla qubits is 
modelled with an additional channel before the encoding op- 
eration. 

as the result of a CPTP channel acting on the state |0) 
simplifies the calculation of channel fidelities for the error 
correction process. 



C. Results 
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The code which maximizes Fq is presented in Figure 
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FIG. 5. Optimal 3-qubit error correcting code against de- 
phasing, when the ancilla qubits used are in mixed states. 



A simple argument illustrates the function of the added 
Toffoli gate in this code. Assuming that the message 
qubit begins in the pure state l^), and that no error oc- 
curs while the encoded state is being transmitted. With- 
out the added Toffoli gate, the unmodified three-bit code 
falsely detects an error on the message qubit and flips 
its state to X \ip). With the added Toffoli gate, the mes- 
sage qubit is flipped before it is encoded, and the error- 
correcting gate at the end restores the original state 

The fidelities of both the traditional three-bit code and 
the optimal code can be expressed in terms of q, the 
mixing parameter introduced in Equation |11| 
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If both ancilla qubits have a polarization 0.6 (q = 1 — e = 
0.4, corresponding to the electron polarization at approx- 
imately 2 K), the difference between these fidelities is 
readily apparent; it is seen in Figure [6] 

Fidelities for these codes can be expressed as functions 
of the ancilla density matrix elements. Thermalization 
is assumed to have occurred before the error correction 
procedure begins, producing a diagonal ancilla density 
matrix: 



Poo|00) (00| + poi |01> (01| 
+ Pxo ]10> <10| + pii |11> <11| 



(12) 



for which F c(Traditional) = pn(2p - 1) + 

(1 -p)(l + (p - 2p 2 )(l - 2( Pl0 + p 01 ))), (13) 

^(Optimal) = (1 - P)(l + (P ~ 2p 2 )(2p Q0 - 1)). (14) 



^(Optimal) is greater than 1-p whenever p Q0 > 1/2. 
This required poo is contrasted with the critical density 
matrix element for correction using the traditional code 
in Figure [7] poo > V 2 w iH serve as the criterion for 
useful error correction. However, the equilibrium room- 
temperature polarization of carbon-13 (a commonly-used 
nucleus) in a 7 T field results in p 00 — l /i + 4.8 x 10~ 7 . 
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FIG. 6. The channel fidelity of three processes: transmis- 
sion through a bit-flip channel, error correction with the tra- 
ditional code and error correction with the code discussed 
above. Unlike the traditional code acting on clean qubits, the 
leading-order error term remains linear, although the q 2 term 
is eliminated. 



Critical p Q0 for Error Correction 




0.2 



— Optimal Code 
Traditional Code 



0.1 



0.2 



0.3 



0.4 



FIG. 7. Ancilla density matrix element required in order for 
error correction to improve the fidelity of a state exposed to a 
noisy channel. This is the criterion which determines whether 
or not the optimal error correction code produces a fidelity 
which is greater than 1 — p, the fidelity produced by inaction. 
Since the fidelity produced by the traditional code does not 
depend exclusively on poo, the matrix element is calculated 
from q, the mixing parameter. 



If refrigeration were used to enhance the polarization 
of these qubits, the temperature required would be 0.6 
mK. It is impossible to adequately polarize these qubits 
with current (or near-future) technology by increasing 
the magnetic field, since 



e = tanh 



\k B T 



TL-i. oc Bo 



(15) 



Doubling the magnetic field present (a technical chal- 
lenge) raises the temperature for constant polarization 



to 1.2 mK. A different technique for enhancing the po- 
larization of a quantum register is, therefore, necessary. 



IV. HEAT BATH ALGORITHMIC COOLING 

A. Background 

Given the need for a different cooling technique, Schul- 
man and Vazirani |15) outlined a procedure to extract a 
small number of highly-polarized qubits from a large en- 
semble with low polarization. Using the closed-system al- 
gorithm presented there, it is possible to extract O (e 2 n) 
qubits with a constant bias from a register of n qubits 
with initial polarization e. To extract the required three 
qubits with polarization v2— 1 (which is required in order 
to satisfy the requirement that poo > V 2 ) from a regis- 
ter with an initial polarization of ~ 1CU 6 would require 
a molecule with ~ 6 x 10 13 spins, clearly an impractical 
criterion. 

Open system methods provide a much better prospect 
for implementation. An efficient method has already 
been implemented at the date of this writing, the partner- 
pairing algorithm introduced by Schulman et al |16j and 
implemented by Ryan et al. [7]. 



B. The Partner Pairing Algorithm 

The partner pairing algorithm was first discussed in 
[16] . where it was found to be optimally efficient (it pro- 
vides the maximum polarization increase possible per op- 
eration) and scalable (the number of elementary opera- 
tions scales polynomially with the system size). There 
are two elementary operations in the partner pairing al- 
gorithm; compression and exchange. Compression ap- 
plies a unitary operator to the density matrix of a 
thermally-mixed state, arranging the density matrix ele- 
ments so that the element corresponding to |0)® n (0|® n 
is the largest element, and the rest are non-increasing. 
Exchange replaces the last qubit in the thermally-mixed 
register with one from the heat bath; it is the net effect of 
allowing that qubit to come to thermal equilibrium with 
the bath. 

These two operations, combined and iterated, will bias 
the density matrix of the ancilla toward |0)® n . Since the 
states involved are classical, and the operations are sim- 
ple, the effect of this algorithm is calculated numerically 
in matlab. There are two situations in which this tech- 
nique will be used; initial state preparation and ancilla 
refreshing. It is important not to cool the qubit which 
stores the information to be preserved by error correction 
while refreshing the ancilla. Also, note that the partner 
pairing algorithm requires three qubits in order to ob- 
tain polarization beyond that of the bath, because once 
two qubits are polarized to the same ebath, their density 
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matrix is of the form 
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Here, the diagonal density matrix elements are already 
ordered, so there is no way to perform a compression 
step. The temperature required for a bath polarization 
needed to perform useful error correction (poo > V 2 ) ^ s 
then 3.4 Kelvin, lower than what can be achieved by im- 
mersing the NMR sample in liquid helium. Therefore, 
it is necessary for experimental feasibility to introduce a 
fourth qubit to be used to refresh ancilla qubits. This 
qubit can also be used to enhance the results of state 
preparation. Since the compression step of each HBAC 
iteration leaves the qubits in the register being cooled 
at different polarizations, we choose to use the most po- 
larized qubits in the error-correction ancilla, and reserve 
the most polarized remaining qubit for the message to be 
transmitted. The purity of this qubit is shown in Figure 

EH 




FIG. 8. The effect of the partner pairing algorithm on the di- 
agonal elements of an initially highly mixed thermal state. 
The even-numbered operations are exchanges between the 
register and the polarization bath, and the odd-numbered op- 
erations are unitary operations which sort the density matrix 
elements according to size. 



C. Results 

In this and the following section, we describe the physi- 
cal parameters required to perform multiple rounds of er- 
ror correction in an experiment and assess the feasibility 
of such an experiment. 

The most important use of algorithmic cooling in this 
proposal is to refresh ancilla qubits between rounds of 
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FIG. 9. The ancilla density matrix element after refreshing 
between error correction steps. A value above 0.5 indicates 
that useful (Fc > 1 — p) error correction can take place. 



error correction. (Algorithmic cooling will also be used 
to prepare the initial states used in error correction, but 
refreshing is the critical step for performing multiple it- 
erations.) In order to obtain the poo required for error 
correction, we subject three qubits to the partner-pairing 
algorithm, and use the two most polarized for the ancilla, 
assuming that the remaining qubit will be thermalized 
before the next refreshing stage. (This qubit is still nec- 
essary to facilitate polarization of the other two, as de- 
scribed in the previous section). We present in Figure [9] 
the achievable poo as a function of the temperature and 
number of iterations of the partner-pairing algorithm to 
be performed on the register. From this, we see that 
temperatures above 4.7 K will not produce usable ancilla 
qubits. Any temperature below this will be useful, with 
colder temperatures being preferable. Also, we see that 
there is very little beneficial effect to performing more 
than three iterations of algorithmic cooling, as long as 
the temperature is below 4.7 K. Since the aim of algo- 
rithmic cooling is to produce polarization on two qubits, 
initialization with the same number of iterations and the 
same temperature will result in a third qubit with con- 
siderable polarization. This polarization is displayed in 
a similar contour map in Figure [TO) This polarization is 
uncommonly large for NMR experiments, increasing the 
efficacy of the protocol as compared to a typical NMR 
experiment. 

It is important to note that the number of qubits which 
are readily available for experiments in NMR is limited. 
With this in mind, two algorithms are proposed; one 
which contains all of the necessary stages to execute mul- 
tiple rounds of error correction with only four qubits, 
and one which has six qubits, so that the message state 
need not be unencoded while ancilla qubits are being 
refreshed. The four-qubit code recommended for exper- 
imental implementation is presented in Figure |11| The 
six-qubit code recommended for practical extraction of 
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FIG. 10. The density matrix element corresponding to |0) (0| 
on the message qubit immediately following initialization. 



entropy is presented in Figure [T2| We determine that, 
given perfect control of the unitary operations we im- 
plement, that three iterations of HBAC will be required 
(the number implemented in [7]) at a temperature below 
4.7 K, which can be achieved in NMR. Coherent control 
over electron degrees of freedom will also be required, 
which is experimentally feasible. In the next section, we 
re-assess the feasibility of an NMR experiment to per- 
form multiple rounds of error correction, given the fact 
that implementation of unitary operators in experiment 
introduces quantifiable noise. 



D. Effects of Imperfect Control 

Since the goal of error correction is to reduce the effects 
of decoherence, it is important to examine the decoherent 
effects of imperfect implementation of the unitary oper- 
ations used for error correction. In [T7], it is found that 
the average error per gate implemented on malonic acid 
is 1%. This figure can be used to further determine 
the utility of the optimal error-correcting code, and the 
exterior temperature required for continual iteration of 
the error-correction and refreshing procedures. 

First, we examine the effect of imperfect control on 
the error correction unitaries; encoding and decoding the 
fault-tolerant state. In order to model imperfect uni- 
taries, we create a channel consisting of the desired uni- 
tary, followed by depolarization [IB] on all three qubits 
with error parameter c. 



U = 




(16) 



We represent the efficacy of a given quantum operation 
using the channel fidelity of the depolarizing map on 
three qubits: 
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We plot in Figure [13] the density matrix element required 
to perform useful error correction at a series of estimated 
gate fidelities. 

We now examine the impact of imperfect control on 
algorithmic cooling itself. We implement a depolarizing 
map after each compression operation, similar to the pro- 
cess presented in Equation [16] In Figure [l4j we present 
the results of refreshing with imperfect control for 6 and 
8 iterations of HBAC, with a bath polarization of 0.31 
and 0.36 (corresponding to temperatures of 4.7 and 4 K, 
respectively). In Figure 15 the same analysis is under- 
taken for initialization. 



V. CONCLUSION 

Here, we have analysed the use of error correction for 
data storage and processing. The transport of entropy 
from a qubit to be preserved to the exterior environment 
is proposed, as is the use of an intermediate heat bath. 
The procedures outlined here for state initialization and 
ancilla re-polarization are experimentally achievable with 
existing ESR/NMR equipment, as well as in other poten- 
tial quantum computing architectures where a heat bath 
is present. 
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FIG. 11. Proof of concept /testbed algorithm for multiple rounds of error correction. After the partner-pairing algorithm is 
used to initialize the register, an arbitrary 1-qubit state is prepared on the qubit with the third-highest polarization, the two 
most polarized qubits being used as an ancilla. Within the dotted rectangles, fault-tolerant computations can be implemented. 
Note that the message qubit is left unencoded while the ancilla is being refreshed. The fourth qubit is employed only in the 
partner-pairing algorithm in order to obtain polarization in excess of the bath polarization. Since the number of iterations in 
the partner-pairing algorithm is variable, its instances in the algorithm above have been left as subroutines. 
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FIG. 12. Plausible algorithm for storage and computation on one logical qubit. After the partner-pairing algorithm is used to 
initialize the register, an arbitrary 1-qubit state is prepared. Here, an extra two-qubit register is used for error correction while 
the other 2-qubit ancilla is being refreshed. If, in the experiment chosen, the partner pairing algorithm can be implemented 
on non-adjacent qubits (as is the case in solid state NMR), then the two swap gates between decoding and encoding can be 
eliminated. The message qubit can then be encoded at all times, and the chief source of error is imperfect control. 
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FIG. 13. Ancilla density matrix element required to perform 
useful error correction, analogous to Figure [7| Curves based 
on gate fidelities given in the legend are shown. Note the 
presence of the curve with 99% gate fidelity, this results from 
the estimated 1% error rate for single-qubit gates in NMR. 
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FIG. 14. Polarization produced by 6 iterations of imperfect 
HBAC on 3 qubits, with an exterior bath polarization of 0.3. 
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FIG. 15. Polarization produced by 6 iterations of imperfect 
HBAC on 4 qubits, with an exterior bath polarization of 0.3. 



